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Abstract 

The design process associated with large engineering 
systems requires an initial decomposition of the 
complex system into modules of design tasks which are 
coupled through the transference of output data, hi 
analyzing or optimizing such a coupled system, it is 
essential to be able to determine which interactions 
figure prominently enough to significantly affect the 
accuracy of the system solution. Many decomposition 
approaches assume the capability is available to 
determine what design tasks and interactions exist and 
what order of execution will be imposed during the 
analysis process. Unfortunately, this is often a 
complex problem and beyond the capabilities of a 
human design manager. A new feature for DeMAID 
(Design Manager's Aid for Intelligent Decomposition) 
will allow the design manager to use coupling strength 
information to find a proper sequence for ordering the 
design tasks. In addition, these coupling strengths aid 
in deciding if certain tasks or couplings could be 
removed (or temporarily suspended) from consideration 
to achieve computational savings without a significant 
loss of system accuracy. New rules are presented and 
two small test cases are used to show the effects of 
using coupling strengths in this manner. 

Introduction 

Many design projects are large and multidisciplinary in 
nature. Before the design of complex systems can 
begin, a significant amount of time and money must be 
invested in determining the order of and the interactions 
among the design tasks. This is particularly true if each 
task is performed by a separate group of engineers or if 
the multidisciplinary analysis is to be incorporated 
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within an optimization procedure. In 1989, a new 
knowledge-based tool [1] was developed to aid the 
design manager in ordering these tasks, determining 
their interactions, grouping iterative processes, and 
displaying the entire process in the format of a design 
structure matrix (NxN). This tool is called the Design 
Manager's Aid for Intelligent Decomposition 
(DeMAID). Since its release to the public, new 
enhancements have been incorporated into DeMAID [2], 
Now that DeMAID has been distributed to several 
industries and universities, feedback from surveys 
indicate that DeMAID users are not only interested in 
just ordering the design process, but want to use 
DeMAID as an aid in reengineering the design problem 
around the flow of information and to emphasize the 
tracking of cost and time. 

Much of the cost and time of a design project is incurred 
in the iterative processes called circuits. These are tight 
groupings of tasks which must be repeated until the 
outputs converge to stable values. DeMAID currently 
groups all the tasks in their respective circuits and 
orders the tasks within the circuit to minimize the 
number of feedback loops. This ordering may or may 
not be the most efficient in terms of cost and time, 
since several different orderings within a circuit may 
yield the same number of feedbacks. One method for 
making the ordering more efficient is to take advantage 
of the strengths of the couplings (interfaces) among the 
tasks [3]. This paper gives a brief overview of new 
directions for DeMAID, followed by a discussion of 
ordering the design process based on coupling strengths. 

New Directions for DeMAID 

When DeMAID was first developed, the knowledge base 
was written around rules for organizing 
multidisciplinary design and optimization problems 
with design variables, behavior variables, constraints, 
and an objective function. As more people started to 
use DeMAID, it quickly became apparent that the tool 
was applicable to a much broader range of problems 
requiring new capabilities. In addition to just ordering 



the design tasks, determining the interfaces, grouping 
the tasks into iterative processes called circuits, and 
decomposing the circuits into a hierarchical display, 
design managers have a need to see what tasks are 
affected if a change is made to some data or task in the 
design process. In addition to tracking the flow of 
tasks, they want to be able to track the flow of 
information through the design process. They want to 
track costs and time in design, manufacturing, and 
maintenance; and use DeMAID to aid in determining 
where simplifications might be made. 

To aid in determining the effects of a change in the 
design process, new rules for selecting tasks for re- 
execution were added to DeMAID [2] . As an example. 
Figure 1 displays a portion of a design structure matrix. 
In this figure, the boxes represent modules, a horizontal 
line from a box represents an output produced by that 
module, and a vertical line to a box represents an input 
required by that module. A dot indicates a coupling 
between two modules because the output of one module 
(e.g. module 25) provides input to other modules (e.g. 
modules 27 and 28). If a change is made to the input of 
module 5, a design manager may want to know what 
effects it would have on the output of module 29. 
Would all modules between module 5 and 29 have to be 
re-executed or only a portion? 



Figure 1. Design structure matrix. 

Figure 2 shows that only a portion of the modules 
(those that are not shaded) require re-execution. 
Although the output from module 18 is an input to 
module 29, this output is not affected by a change to 
module 5. Modules 23, 24, and 26 may be affected by 
the changes to module 5, but since their output does not 
change any input to module 29, they too may be 
omitted from re-execution. Since one or more of these 
tasks may require substantial resources in cost, time. 


personnel, and/or equipment, any task not having to be 
re-executed may result in a significant savings. 



Figure 2. Modified design structure matrix. 

To enable the design manager to track more detailed 
information requires a change in the way the output is 
represented in the format of a module definition. 
Initially [1], the format was: 

(module number name type time output 
status input-list) 

where the output field signifies all output produced by a 
module, whether it be a single piece of data, a string of 
data, a vector, a matrix, or any combination. Design 
managers have expressed their desire to see more details 
than the present representation allows. 

For example, suppose a module represents a finite 
element analysis program which computes stresses, 
displacements, vibration modes, etc. In the original 
format, all this output would be lumped into a single 
name (e.g. feout ). 

(module 1 ml 1 25 feout uk fein geom) 

In Figure 3, this task is represented as module 1, and 
has interfaces to modules 2, 3, and 4. Hence, these 
modules require feout as input. Suppose however that 
module 2 really only requires stresses as input to 
compute stress constraints, module 3 only requires 
displacements to compute displacement constraints, and 
module 4 requires stresses, stress constraints, and 
displacements. The design structure matrix would 
appear the same, but if a change occurred that only 
affected the stresses, then module 3 would not have to 
be re-executed since only the displacement data from 
feout are used as input. 
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(module 1 ml 1 25 feout uk) 

(module 2 m2 1 10 strcons uk) 

(output feout stress disp) 

(requires uk VS strcons stress) 

Figure 3. Modules with multiple outputs. 

The design manager can now trace specific pieces of data 
and be as detailed as the problem requires. This is 
accomplished by adding two new sets of lists (examples 
in Figure 3) to the DeMAID input. The user creates the 
module list as before, but omits the input-list. 
DeMAID uses the two new lists to create the input-list 
for each module. 

The first new list is an output list in the format: 

(output name output-breakdown) 

where the name in the output list matches the output in 
the module list. The output-breakdown in the output 
list is a breakdown of all the pieces of output data the 
design manager needs to trace. 

The second new list aids in determining the interfaces. 
It has the format: 

(requires status strength output input) 

where the output and input fields are elements of 
output-breakdown from the output lists and the strength 
field represents the coupling strength between those two 
elements. 

The concept of coupling strengths becomes very 
important when trying to reduce costs and cycle time. 
Since much of the cost and time of a design project are 


incurred in the iterative processes called circuits, these 
coupling strengths can aid in deciding if certain tasks or 
couplings can be removed (or temporarily suspended) 
from consideration to achieve computational savings 
without a substantial loss of system accuracy. 
DeMAID handles seven levels of coupling strengths: (1) 
extremely weak (EW), (2) very weak (VW), (3) weak 
(W), (4) nominal (N), (5) strong (S), (6) very strong 
(VS), and (7) extremely strong (ES). Currently, the 
user must estimate and input these coupling strengths if 
the requires lists are needed. The sections below 
describe a method for computing the coupling strengths 
as well as rules for using them to reduce time and costs 
with minimal loss in accuracy. 

Coupling Strengths 

Coupling Strength Definition 

Here, coupling strengths are defined in terms of the 
local normalized sensitivities. Since these local 
sensitivities are used in the Global Sensitivity 
Equations (GSE) to obtain total behavioral response 
derivatives, they are already available to the design 
manager [4], 

The GSE approach involves using the chain rule to 
define the total derivatives of the output response 
quantities in terms of local sensitivities of each 
subsystem. The GSE's feasibility has been 
demonstrated in previous work [4,5,6], Once the local 
sensitivities are known, the total derivatives of the 
output response quantities with respect to the design 
variables can be determined from the solution of the 
matrix set of global sensitivity equations. Since the 
dY/dX information is invariably required in design for 
trade-off studies, as well as formal gradient-based 
optimization, it is inevitably available during the design 
process. 

Since the components of the output response vector Y 
and the design variable vector X are of varying 
magnitudes from one subsystem to another, it is 
necessary to scale the local sensitivities. Hence, a 
normalization scheme [7] can be implemented to ensure 
that the conditioning of the system is such that accuracy 
of the solution is not threatened. The local normalized 
sensitivities (denoted by ') for a sample two subsystem 
(A and B) problem are: 

dYA YB dYA , dYB YA dYB 

= and = 

dYB YA dYB dYA YB dYA 
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Sensitivity-Based Heuristics 


The local subsystem sensitivity information can be used 
to quantify the strengths of participating analysis 
couplings. This coupling information can then be used 
to provide the basis for developing heuristics that 
indicate which couplings are “weak” enough to be 
temporarily or permanently suspended. For example, if 
it is determined, from the system in Figure 4a, that the 
coupling from module 4 to module 3, and from module 
5 to module 3 are negligible (i.e. dY3/dY4' and 
dY 3/ dY5' are extremely small in comparison to other 
normalized local sensitivities), it may be decided to 
temporarily suspend the feedback loop among these 
modules, thereby reducing the complexity of the 
problem to that seen in Figure 4b. 



Figures 4a and b. Subsystem modules with full and 
modified dependencies. 

Obviously, in a complex problem involving 
computationally expensive analyses (such as structural 
finite element analyses) that must be executed within an 
iterative framework, the ability to reduce the system 
complexity without sacrificing solution accuracy is of 
the utmost importance. 

The question then becomes, “to what extent may 
solution accuracy be compromised in order to achieve 
solution efficiency?” It is this question that is addressed 
in the development of new rules for DeMAID. The 
application of a knowledge -based system enables the 
answer to this question to be largely problem dependent, 
thus increasing applicability and flexibility. 

Coupling Strength Quantification 

In this paper, a linear distribution between upper and 
lower bounds of the local sensitivity space is used to 
quantify the seven levels of coupling strengths. 

One possible approach is to calculate the mean and 
standard deviation where the mean value of the local 
normalized sensitivity derivatives (sj) can be determined 
from 


s 


1 

N 


N 


E 


i=l 


where N is the number of couplings. The associated 
standard deviation can now be determined from the 
relation. 


c(s) 


— - — L(sj 
N - 1 ;=i V 



— il/2 


The upper and lower bounds of the local normalized 
sensitivity space are defined in terms of the mean value 
and standard deviation as. 


s u = s + k x o( s) 
s 1 = s - k 2 <J( s) 


where l< ( and k 2 are user-prescribed values based on 
experience and heuristics. Anything outside these 
bounds is either extremely weak (EW) or extremely 
strong (ES). 

Once these bounds are defined, associated coupling 
strengths can be assigned to all other couplings based 


on a linear distribution between the bounds. 

Hence, 

if 

Si > s + k,cr(is) 

then 

S( is ES 

if 

3 

s + — k,cr(s) < Sj < s + kjCT(s) 

then 

Sj is VS 

if 

1 3 

s + — k,cr(s) <Si<s + — kjCrii) 

then 

Sj is S 

if 

s - ^-k 2 cr(s) <Sj < s + -^kjOfs) 

then 

Sj is N 

if 

3 1 

s -—k 2 cr(s) <Sj < s ~—k 2 (j(s ) 

then 

Sj is W 

if 

3 

s - k 2 o{s ) <Sj< s — — k 2 o{s) 

then 

SjisVW 

if 

Sj < s — k 2 o(s) 

then 

Sj is EW 


Another possible approach to quantifying these 
coupling strengths is more heuristics-based. Once the 
mean is determined as above, an upper and lower limit 
on the normalized sensitivity space can be prescribed by 
the design manager. The range for the nominal 
quantification can also be pre-determined, leaving a 
linear distribution for the quantifications between the 
upper limit and the nominal range and those between 
the lower limit and the nominal range. For example, if 
the upper bound on the normalized sensitivity space is 
UB, the lower bound is LB, and the distance from the 
mean to the nominal boundary is NB, the 
quantifications would be classified as. 
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if s, > UB 

( s + NB) + UB ^ ^ TTD 
if < s t < UB 


if s + NB < Si < 


(s + NB) + UB 


if s - NB < Sj <s + NB 

, (s — NB)- LB 
if < Sj < s - NB 


if LB < Sj < 
if Sj < LB 


(s - NB) — LB 


then 

Sj is ES 

then 

Sj is VS 

then 

Sj is S 

then 

Sj is N 

then 

Sj is W 

then 

Sj is VW 

then 

Sj is EW 


The coupling strength quantifications are then supplied 
to DeMAID in the requires list. 


modules believed to be important for the complex 
system. DeMAID would yield an ordering of those 
participating tasks based on its planning and scheduling 
components. A full analysis would then be performed 
followed by determination of sensitivity information. 
From the local normalized sensitivities, coupling 
strengths would be quantified and categorized as 
extremely weak (EW), very weak (VW), etc. They 
would then be supplied to DeMAID to determine the 
reduced system and its associated ordering of tasks. The 
reduced system and it's optimization problem would be 
solved for some number of prescribed cycles, n, after 
which the full system and new coupling strengths 
would be evaluated. This approach results in the 
potential for substantially reduced computational times 
with minimal loss in system and optimization accuracy. 


DeMAID Coupling Strength Rules 

Once the coupling strengths have been determined, they 
can be used with a new set of DeMAID rules to 
determine which modules may possibly be removed (or 
temporarily suspended) from the design process. The 
rules work in the following manner. First, all modules 
with at least one coupling having a nominal (N) or 
stronger strength are retained. Modules with only 
extremely weak (EW) coupling strengths are removed. 
The remaining modules are checked to determine their 
relationship to other modules in the system. The 
maximum coupling strength is found for each module 
remaining in the system. If a module has a weak (W) 
maximum coupling strength and interfaces with a 
module with a very strong or extremely strong (VS or 
ES) coupling strength, then that module is retained. 
Otherwise, it is removed from consideration. If a 
module has a very weak (VW) maximum coupling 
strength and interfaces with a module with an extremely 
strong (ES) coupling strength, then that module is 
retained. Otherwise, it too is removed from 
consideration. All references to the output generated by 
the removed modules are deleted from the input lists of 
the remaining modules. The user can now input this 
reduced set of modules into DeMAID. 

A pplication to Analytical Systems 



Figure 5. Interactions in two subsystem non-hierarchic 
environment. 


Intelligent Decomposition 

Before describing the application of the new DeMAID 
rules to analytical systems, it is necessary to understand 
how coupling strength information would be used in a 
design procedure. Figure 5 demonstrates how DeMAID, 
together with the coupling strength data, can 
significantly reduce the number of both full system 
analyses as well as full optimization runs. The design 
manager's judgment would be used to identify all 


Two analytical problems have been developed to 
emulate a complex engineering system with lateral 
couplings. Figure 6 demonstrates the coupling for the 
first analytical system. Within each subsystem, 
outputs are associated with contributing functions. For 
example, within subsystem 1 of Case 1 (Figure 6), 
there are two functions associated with Subsystem 1, 
S1W1 with output wl and S1W2, with output w2. 
These outputs, which may be individual or multiple 
pieces of data, then form input to other subsystems, and 
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hence, functions within those subsystems. The same is 
true for the system of Case 2. 


wl,wi 



zl,z2,z3 


Subsystem 2 
(Output - y) 



yl,y3 wl,w2^ 

yl,y2,y3,y4 
zl,z2,z3 


Figure 6. Non-hierarchic analytical system providing 
testbed for coupling strength comparisons. 

Analytical Coupled System Analyses 


S 1X2 
S2X 1 
S2X2 
S2X3 
S3X1 
S3X2 
S2V3 
S1N1 
S3Z2 
S 1N2 
S2V 1 
S2V4 
S2V2 
S3Z1 
S3Z3 



Implicit relations are defined for both cases for each 
subsystem. The notation f(. . .) denotes a functional 
relationship among the listed inputs and outputs. The 
explicit equations for Case 1 are as defined in [3], while 
those for Case 2 are listed in the Appendix. 

Case 1: 

Subsystem 1 (SSI) 


XI = ( xll ,xl2) and 

wl =f( XI, y3 , zl , z3 ) 

w2 =f( XI , yl , y3 , z2) 

Subsystem 2 (SS2 ) 

X2 

= ( x21 , x22 , x23 ) and 

yi 

= f( X2 , wl , w2 , z2 ) 

y2 

= f( X2 , w2 , zl , z2 ) 

y3 

= f( X2 , zl , z3 ) 

y4 

= f( X2 , wl , z2 ) 

Subsystem 3 (SS3 ) 

X3 

= ( x31, x32 ) and 

zl 

= f( X3 , w2 , y2 , y4 ) 

z2 

= f( X3 , wl , yl ) 

z3 

= f( X3 , y3 , w2 , y4 ) 


The application of DeMAID to the above system of 
coupled functions results in the ordering shown in 
Figure 7. Each function is denoted by a reference to its 
subsystem and to the output associated with it. 
Modules pertaining to the design variables (e.g. S1X1) 
are included to make the identification of subsystem 
inputs easier. 


Figure 7. Case 1 modules with minimized feedbacks. 


Case 2: 


1 (SSI ) 


XI 

= ( xll , xl2 , xl3 ) and 

tl 

= f( XI, yl , v2 ) 

t2 

= f(Xl , vl , zl ) 

t3 

= f( XI , v2 , z2 , u2 ) 

t4 

= f( XI , w2 , v4 , y3 ) 

Subsystem 2 (SS2 ) 

X2 

= ( x21 , x22 ) and 

ul 

= f(X2 , zl , w3 ) 

u2 

= f( X2 , tl , t3 , w2 , z3 ; 

u3 

= f( X2 , tl , yl ) 

Subsystem 3 (SS3 ) 

X3 

= ( x31, x32 , x33 ) and 

vl 

= f( X3 , t2 , u2 , zl , w2 

v2 

= f( X3 , tl , y2 ) 

v3 

= f( X3 , u3 , wl ) 

v4 

= f( X3 , yl , y3 , t4 , w3 

Subsystem 4 (SS4) 

X4 

= ( x41 , x42 , x43 ) and 

wl 

= f( X4 , yl , y2 , z2 ) 

w2 

= f( X4 , zl , z2 , z3 , u2 

w3 

> 

”2; 

X 

n 

Subsystem 5 (SS5 ) 

X5 

= ( x51, x52 , x53 ) and 

yi 

= f( X5 , tl , u3 , v3 ) 

y2 

= f( X5 , u3 , wl , v2 ) 

y3 

= f( X5 , z3 , t4 , ul ) 

y4 

= f( X5 , t3 , u2 , z3 ) 

Subsystem 6 ( SS6 ) 

X6 

= ( x61 , x62 ) and 

zl 

= f( X6 , t2 , t3 , vl ) 

z2 

= f( X6 , v2 , v3 , wl ) 

z3 

= f( X6 , wl , w2 , 


t2 , u2 


y4 ) 
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Figure 8 shows the full system for Case 2, with all 
couplings. 



Figure 8. Case 2 modules with minimized feedbacks. 

Discussion of Results 

The first step in applying the DeMAID-based reduction 
scheme is to quantify the coupling strengths for Cases 1 
and 2. The results of applying the heuristics-based 
quantification scheme previously described are shown in 
Table 1 for Case 1 and Table 2 for Case 2. The 
distribution for the normalized local sensitivity space is 
shown in Figures 9 and 10 for the two cases. 


.0002 .0647 .1292 .6292 .8146 1.00 



EW VW W N S VS ES 


Figure 9. Case 1 boundaries for quantifying coupling 
strengths. 


.00001 .04901 .0980 .4980 1.249 2.00 



EW VW W N S VS ES 


Figure 10. Case 2 boundaries for quantifying coupling 
strengths. 

For Case 1, the new rules resulted in modules S3Z2 and 
S3Z3 being removed from consideration as in the 
original paper presenting this case [3], As reported in 
that paper, the largest percent difference, 7.153%, was 
in the solution for w2. The next largest, 2.28% was 
associated with yl, and all others were below 1%. 


For Case 2, the same set of rules allowed modules v3, 
w2, y4, and z3 to be removed. As can be seen from 
Table 3 the output values for the reduced system are 
very close to those of the full system. The removed 
modules are shaded in the table. Removing the 
modules, however, produced an interesting change in the 
design structure matrix as shown in Figure 11. One 
module was removed from the top circuit, reducing the 
feedbacks by one. The bottom circuit was unchanged. 
Three modules were removed from the middle circuit 
resulting in a division into two smaller circuits. 



Figure 1 1 . Case 2 after modules removed 

With these two cases, it has been shown that general 
rules can be developed to remove (or temporarily 
suspend) modules from execution by examining the 
strengths of the couplings among the modules. The 
strengths or weaknesses of individual couplings can be 
determined by applying statistical techniques to 
sensitivity data. With this knowledge, modules may be 
removed from the design process which offers a 
potential improvement in both time and cost without a 
loss in solution accuracy. 
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Table 1. Coupling Strength Quantifications for Case 1 using heuristic approach. 



delwl/ 

delw2/ 

delyl/ 

dely2/ 

dely3/ 

dely4/ 

delzi/ 

delz2/ 

delz3/ 

delwl 



ES 



W 


VW 


delw2 



N 

VW 



VW 


W 

delyl 


VW 






vw 


dely2 







VW 



dely3 

VW 

VW 







VS 

dely4 







s 


EW 

delzi 

vw 



VW 

VW 





delz2 


w 

VW 

EW 


EW 




delz3 

EW 




VW 






Table 2. Coupling Strength Quantifications for Case 2 using heuristic approach. 
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Table 3. Comparison of output values for 


full and reduced systems for Case 2. 


Output 

Full System 

Reduced System 

1 1 

71.3021 

71.3019 

t 2 

31.5566 

31.5560 

1 3 | 

94.8302 

94.8267 

t 4 

38.7446 

38.7472 

ul 

70.3609 

70.3600 

u2 

90.1629 

90.1616 

u3 

37.4936 

37.4936 

V 1 

135.3255 

135.3257 

v 2 

95.4775 

95.4776 

v 3 

143.4225 

143.4225 

v 4 

119.3493 

119.3497 

w 1 

141.7335 

141.7328 

w 2 

34.4332 

34.4332 

w 3 

132.2914 

132.2894 

y i 

139.2781 

139.2783 

y 2 

81.6617 

81.6616 

y 3 

147.3626 

147.3626 

y 4 

95.5950 

95.5950 

Z 1 

141.3610 

141.3611 

z 2 

77.4906 

77.4906 

z 3 

110.3259 

110.3259 


Concluding Remarks 


Much of the cost and time of a design project are 
incurred in the iterative processes called circuits. A 
new method for making the ordering within the circuits 
more efficient by taking advantage of the strengths of 
the couplings among the tasks was introduced. This 
new method applies new rules in the knowledge base 
based on coupling strength information derived from 
sensitivity data. These rules determine which modules 
may be removed (or temporarily suspended) from the 
design process without sacrificing solution accuracy. 
Two coupled systems composed of analytical equations 
were tested for verification purposes using the same 
rules. In both cases, modules were removed based on 
the sensitivity data with little or no effect on the 
solution. 
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Appendix 


Explicit analytical expressions for each subsystem analysis of Case 2 are presented. Each output quantity (e.g. wl) 
is assumed to represent output from a subroutine. Table A1 shows coefficient values used in the system analysis. 


Subsystem 1 (SSI) 


XI 
tl 
t2 

t3 
t4 

Subsystem 


X2 
ul 
u2 
u3 

Subsystem 


(xll , xl2 
f( XI, yl , 
f( XI , vl , 
f( XI , v2 . 
f( XI , w2 
2 < SS2 ) 


X3 
vl 
v2 
v3 
v4 

Subsystem 


( x21 , x22 
f( X2 , zl , 
f( X2 , tl , 
f( X2 , tl , 
3 (SS3) 


, xl3 ) and 
v2 ) 
zl ) 

, z2 , u2 ) 

, v4 , y3 ) 

) and 
w3 ) 

t3 , w2 , z3 ) 

yi ) 


= al(xll) 2 + a2(xl2) 1/2 + a3(yl) 1/2 - a4(v2) 

= a5(xl2) + a6/xl3 - a7(vl) + a&zl) 2 

= a9(xl 1) -fc,aJ0/xl2 + al l(xl3) 2 - al2(v2) 2 - al3(z2) 2 - al4(u2) 

= al5(xl 1) 3/2 + al6(xl2) + al7(xl3) 2 + al8(w2) 172 - al9(v4) + a20/y3 


bl(x21) z + b2(x22) + b3(zJ) - b4(w3) 1/2 
- b5(x21) + b6/x22 - b™ti vZ a. ^ h c 

= bl l(x22) 2 + bl2(tl) 


b5(x21)_+ 0 b6/x2_2 - b^tlj^+ b8(t^ 0 + b9/(w2) 1/2 - bl0(z3) 1/2 


X4 
wl 
w2 
w3 
Subsystem 


( x31, x32 
f( X3 , t2 , 
f( X3 , tl . 
f( X3 , u3 
f( X3 , yl . 
4 (SS4) 


u2 33 zl fwi )=cl(x31) + c2(x32) 2 - c3(t2) + c4(u2) 1/2 + c5(zl) - c6/(w2) 1/2 
y2 ) = c7(x32) 1/2 + c8(x33) T c9(tl) 2 + cl0(y2) 

, wl) = cl l(x31) 2 + cd2(x33)l /2 + cl3(u3) r/2 - cl4(wl) 

y3 , t4 , w3 , ul )= cl5(x32) z - cl6(yl) + cl7/(t4) + cl8(w3) I/2 + cl9(ul) - c20/(y3) 1/2 


X5 

y i 
y 2 
y 3 

y4 

Subsystem 6 ( SS6 ) 


( x41 , x42 
f( X4 , yl . 
f( X4 , zl , 
f( X4 , t4 , 
5 ( SS5 ) 

‘ ' x52 
, tl , 
, u3 . 
, z3 . 
, t3 , 


(x51. 
f( X5 
f( X5 
f( X5 
f( X5 


X6 

zl 

z2 

z3 


(x61 , x62 
f( X6 , t2 , 
f( X6 , v2 . 
f( X6 , wl 


, x43 ) and 
, y2 , z2 ) 
z2 , z3 , u2 )= 
v4 ) 

, x53 ) and 
u3 , v3 ) 

, wl , v2 ) 

, t4 , ul ) 
u2 , z3 ) 

) and 
t3 , vl ) 

, v3 , wl ) 

, w2 , t2 , u2 


= dl(x41) 3 + d2(x42) + d3(yl) - d4(y2) 2 + d5/z2 
= d6(x42) 2 - d7(zl) - d8(z2) + d97z3 + dl0(u2) 

= dll(x4l) + dl2/x42 + dl3(x43) 2 - dl4(t4) - dl5(v4) 

= el(x51)J- e2(x52) + e3(x53) 2 - e4(tl) 2 + e5(u3) + e6(v3) 1/2 
= e7(x51) 2 + e8(x53) + e9(u3) + elO/wl - ell(v2) 

= el2(x52) + el3/x53 d- el4(z3) - el5/t4 - el6(ul) 

= el7(x51) 2 + el8(t3) 2 - el9(u2) + e20/z3 

= fl(x61).+ f2(x62) 2 - f3(t2)= f4(t3) + f5.jC.vl ) 

= f6(x61) 2 + f7/x62 + f8(v2) 2 - f9(wl) 1/2 - fl0/v3 
y4 )= fl l(x62) - fl2(wl) + fl3(t2) 172 + fl4(u2) + fl5/w2 - fl6(y4) 1/2 


Table Al. Coefficient Values for Analytical System of Case 2. 


Coefficients 

Coefficient values 

al-a20 

23.521 45.122 .9352 .2835 23.175 30.645 .5392 .0024 12.944 43.846 18.345 .0113 .0021 
.1635 12.347 21.478 7.745 .0563 .7256 .9476 

bl-bl3 

8.435 12.563 .1577 .9388 24.547 24.356 1.457 .3846 .9862 .7377 9.367 .0285 .6758 

cl-c20 

16.846 10.360 .3367 .157 .0054 .4783 23.644 14.466 .0064 .5347 8.846 44.467 .3748 .1094 
14.536 .0949 .4503 .1003 .0083 .0776 

dl-dl5 

14.896 13.746 .4567 .0144 .0944 16.436 .0113 .0673 .0356 .0275 8.638 26.693 6.536 .6685 
.8467 

el-e20 

27.783 35.552 10.377 .0059 .0983 .0633 15.367 18.653 .1079 .2763 .831 22.561 35.649 
.6463 .7366 .0667 9.882 .0011 .0359 .2033 

fl-fl6 

8.629 7.932 .0274 .0235 .0021 14.367 48.837 .0009 .0328 .0825 29.663 .0846 .0024 .0436 
.0875 .0289 
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